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Abstract 

o . 

£\j i The Kudryashov-Sinelshchikov equation for describing the pressure waves in liquid with gas bubbles is studied. 

New exact solution of this equation are found. Modification of truncated expansion method is used for obtaining 
exact solution of this equation. 
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1 Introduction 

G ■ 

Recently, Kudryashov and Sinelshchikov [T] introduced the following equation 

u t + juu x + u xxx - e {uu xx ) x - ku x u xx -vu xx - 8(uu x ) x = 0, (1) 
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■ where 7, e, k, v and S are real parameters. Eq. (TjQ) describes pressure waves in the liquid with gas bubbles taking 
0^ , into account the heat transfer and viscosity [T]. We call this equation the Kudryashov-Sinelshchikov equation. 

It is well known that pressure waves in gas-liquid mixture is characterized by the Burgers-Korteweg-de Vries 
(BKdV) equation and the Korteweg-de Vries (KdV) equation [2H5]. The Kudryashov-Sinelshchikov equation is 
generalization of the KdV and the BKdV equation. Indeed, assuming e = K = (5 = 0we have the Burgers- 
Korteweg-de- Vries equation. In the case ofe = ft = A = <5 = 0we get the famous Korteweg-de Vries equation. 
The aim of this work is to find exact solutions of the Kudryashov-Sinelshchikov equation. 

We know some methods for finding exact solution of ordinary partial differential(ODE) equations. Let us 
note some of them: the truncated expansion method [SHS], the simplest equation method [TDJ[TT], an automated 
tanh- function method [T^] , the polygons method [T3] and the Clarkson-Kruskal direct method [2] . 

For finding exact solution of the Kudryashov-Sinelshchikov equation we use the modification of truncated expan- 
sion method that was introduced in [15J . Using truncated expansion method for finding exact solutions of ODEs 
we obtain the overdetermined system of differential equations. As a rule this systems is difficult to solve. The 
modification of this method allows us to transform this system of differential equations to the system of algebraic 
equations. As a result we have essential simplification of solutions construction procedure. 

This paper organized as follows. In the Section 2 we introduce the method applied. Solitary waves solutions of 
Eq. (JTJ in the partial case of v = S = are described in Section 3. Exact solution of Eq. (JTJ in general case are 
discussed in Section 4. 

2 Method applied 

Let us present the modification of the truncated expansion method |15j . We consider the nonlinear partial differential 
equation in the form 

E[u u u x , . ..,x,t] = 0. (2) 

Using traveling wave 

u{x,t) = y{z), z = kx — wt. (3) 
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from Eq. ^ we obtain the ordinary nonlinear differential equation 

L[y,y z ,...,k,w]=0. (4) 

The modification of the truncated expansion method contains from the following steps |15j . 
The first step. Determination of the dominant term with highest order of singularity. To find dominant terms 
we substitute 

y = z~ p , (5) 

into all terms of Eq. (U). Then we should compare degrees of all terms of Eq. Q and choose two or more with 
the highest degree. The maximum value of p is called the pole of Eq. ((4]) and we denote it as N. It should be 
noted that method can be applied when N is integer. If the value N is noninteger one can transform the equation 
studied. 

The second step. We look for exact solution of Eq. (|4]) in the form 

y = a + aiQO) + a 2 {z)Q(z) 2 + ... + a N Q(z) N , (6) 
where Q(z) is the following function 

The third step. We can calculate necessary number of derivatives of function y. It is easy to do using Maple or 
Mathematica package. Using case N = 2 we have some derivatives of function y(z) in the form 

y = a + a\Q + a 2 Q 2 , 

y z = -a 1 Q + (a 1 -2a 2 )Q 2 +2a 2 Q 3 : (8) 
Vzz = oiQ + (4a 2 - 3ai)Q 2 + (2oi - 10a 2 )Q 3 + 6a 2 Q 4 . 

The fourth step. We substitute expressions ©-© in Eq- ©■ Then collect all terms with the same powers of 
function Q(z) and equate this expressions to zero. As a result we obtain algebraic system of equations. Solving this 
system we get the values of unknown parameters. 

This algorithm can be easily generalized to polynomial differential equation of any order. 

3 Exact solutions of the Kudryashov-Sinelshchikov equation in the 

case of v = S = 

Let us find the exact solutions of the Kudryashov-Sinelshchikov equation. Using scale transformation 

x = x', t = t\ u = -u', (9) 

£ 

the Kudryashov-Sinelshchikov equation is written in the form [1] 

u t + auu x + u xxx - (u u xx ) x - fi u x u xx = 0, (10) 

where a = j/s, f3 = n/e (primes are omitted). Taking the traveling wave ansatz ([3]) into account and integrating 
with respect to z from Eq. (fTU]) we have 

k a 
d - lo y + — y 2 + k\y zz -yy zz ^tiy 2 J=0 (11) 

Here C\ is integration constant. 

The pole order of Eq. (fTT|) is N = —-pj^- One can see that at = —3 or = —4 Eq. (TTT1) have the pole of 
second or first order consequently. So we look for solution of Eq. (jlip in the form 

y = a Q + ai Q + a 2 Q 2 (12) 
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Substituting ([12")) into Eq. (fTTj) and taking into account relations ([5]) we obtain the system of algebraic equations 
in the form 

-2k 3 aj (/3 + 3) =0, 
-2 k 3 a 2 (-2 a 2 /3 - 5 a 2 + 4 ai + f3 ax) = 0, 

fc 

-- (8a 2 2 k 2 + 4fc 2 ai 2 - 26 k 2 a 2 a 1 - 12a 2 fc 2 + 12 k 2 a 2 a Q + 
+4 A: 2 /? a 2 . + k 2 /3 a\ - 8 fc 2 /3 a 2 ai - a a 2 2 ) = 0, 

—A: (5 k 2 a 2 ax + 2 fc 2 ai ao + 2 fc 2 /? a 2 ai — 10 k 2 a 2 &o — 2 k 2 a\ — 

-k 2 /3 a{ - 3 fc 2 a 2 + 10 a 2 fc 2 - a ax a 2 ) = 0, ( 13 ) 

kaa\ k 3 f3a\ o , 

— uj a 2 3 fc ax — 4 k a 2 ao + ka ao a 2 + 

+3 fc 3 ai a + 4 k 3 a 2 - k 3 a\ = 0, 
— ax (to — ao a k + ao k 3 — k 3 ) = 0, 

From (fT3")) we have following values of coefficients ao, ax, a 2 and paraments a, 0, uj, C\ 

(fc 3 -w) (2wfc 2 -wci-fc 3 a) 

= — 3 , C\ = o , 

2fc(fc 2 -a) 2 

12fc (fc a — uj ) 12k(ui-ka) k 3 -tu (14) 

a2 ~ k^-o? ' Ql ~ fc 4 - a 2 ' 00 " fc (fc 2 - a) ' 

a ^ ±fc 2 

/3 = — 3, a 2 = -ai a = fc 2 , ao = — + 1, 

6 

,3 n _ fc 3 («l-36) 



uj = k 3 , Ci = 



(15) 



72 

/3 = —3, a 2 = — ax a — — fc 2 , ao = 1, 

„ fc3 (16) 

uj = -fc'\ Ci = 

2 

a 2 = 0, /3 = -4, ai = 2(l-a ), a / 1 

,2 .3 r- a k 3 (2~ao) ( 17 ) 

a = fc , uj = k , Ux = - 



2 

We have four families of solitary wave solutions of Eq. ([TTj) in the form 

k 3 -uj 12k (u-k a) 12fc(fca-w) 

y ~ fc (fc 2 - a) + (fc 4 -a 2 )(l + e z ) + (fc 4 - a 2 )(l + e*) 2 ( ^ 

1 ai , aieZ Cio"\ 

" = 1 - T + (TTPF ( ' 

" = 1+ (IT$P (20) 

' = »» + S iT^ (21) 

corresponding to values of paraments a,/3,w,Ci defined by relations (fT4")l . (fT5")) . (|16l) and (IT7)) . 
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Figure 1: Exact solution of (|18p . solid - k = 1, cj = 0.75, a = 2, dash - k = 1.1, w = 0.75, a = 0.5, dot 
k = 1, u> = 1, a = 0.5. 
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Figure 2: Exact solution of (|T9|). ai = 1.5; 1;0.5 (solid; dash; dot). 



4 Exact solutions of the Kudryashov-Sinelshchikov equation in the 

case of v ^ 0, £ 7^ 

Using transformation (J9j> we can write the Kudryashov-Sinelshchikov equation in the form 

u t + auu x + u xxx - (u u xx ) x - P u x u xx - v u xx - fx(u u x ) x = 0, (22) 

where a = 7/e, /? = k/e, /j. — S/e (primes are omitted). Taking traveling wave transformation ([3]) into account and 
integrating with respect to z from Eq. (|2"21 we have 



Ci - ujy+ '-^-y 2 + k 3 (y zz -yy zz - ^-y 2 z ) - fc 2 (v y z + nyy z ) = 



k a 
~~2 



Let us look for solution of Eq. (f23j) in the form 

y = a + aiQ + a 2 Q 2 



(23) 



(24) 



Substituting (pH)) into Eq. (|2"3")l and with help of relations © we obtain system of algebraic equations from 
which we find values of <Z2, a±, ag, oj, C\, /3 



2{k + v)k A{k + v)k 

p = —6, fJ. = k, a 2 — , a\ = , 

a a 

12k 2 v + 12k 3 -5ak + va a. r , , , 2 ., , 

a = — — , u = -(5 k-v)-k [k + v), 

ba k b 



(25) 



Ci = 



(u-5k) (12 k 2 v + 12 k 3 - 5 a k ■ 
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-3, a = fi 2 -k 2 , 



(l 2 



k 2 (v + n) 
fi(k 2 — fi 2 ) ' 



01 



[y + fi) k 
(k - n)n ' 



a 



(26) 



fc 



2 4fc 2 ) 



a = 



-4, a2 = 0, ai = 



fc 2 + /i fc — a + kv + v fi 



k 2 — a — fi 2 



2(fi + v) k 
[i 2 — k 2 + a 1 

ka (fc 2 — a + v fi\ 
k 2 — a — /i 2 



(27) 



ka (v 2 fi 2 -2aufi- k 2 (p 2 + v 2 ) + (fc 2 - a) 2 



2 (fc 2 - a - fi 2 ) 2 

Using values of parameters (l25t we have following kink-type solution of Eq. (l24l) 

12 fc 2 i/ - 5 fca + 12fc 3 + z/a 4 (fc + v) k 2{k + v)k 
V ~ 6to + a (l+e z ) ~ a (1 + e z ) 2 

With help of relations (|26l) we obtain following solution of Eq. (f24|) 

2 fc - fi + v (fi+u)k k 2 (fi + v) 

V = 2(fc- M ) ~ (fc- M ) M (1+C*) + ^ ( fc 2_ /i 2)( 1+e z)2 

Taking into account values of parameters ([27|) we have following kink-type solution of Eq. ([24 

fc 2 + fik — a + kv + v fi | 2(fi + v)k 



(28) 



(29) 



k 2 -a-fi 2 (p 2 - k 2 + a)(l + e z ) 

Dependence solution (|29j) from z at different values of parameter \i at fc = 2, 2/ = 1 are illustrated on Fig. [3] 

y 



(30) 
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Figure 3: The solution (H| of Eq. pij) at /i = 1; 0.7; 0.4 (curves 1,2,3). 
We believe that solutions (l28l). ([29} and (M are new. 



5 Conclusion 

The Kudryashov-Sinelshchikov equation was studied using the modification of the truncated method. The algorithm 
of the method applied was presented. The efficiency of this method was demonstrated. New exact solution of the 
Kudryashov-Sinelshchikov equation were obtained. 

This work was supported by the federal target programm "Research and scientific-pedagogical personnel of 
innovation in Russia" on 2009-2011, Contracts P 28, P 741. 
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